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^ , Abstract 

in 

The Euler equations of ideal gas dynamics posess a remarkable nonlinear involutional symmetry 
which allows one to factor out an arbitrary uniform expansion or contraction of the system. The 
nature of this symmetry (called by cosmologists the transformation to supercomoving variables) is 
\ discussed and its origin clarified. It is pointed out that this symmetry allows one to map an explosion 

problem to a dual implosion problem and vice versa. The application to laboratory simulations of 
supernova remnants is considered; in principle this duality allows the complete three-dimensional 
evolution of highly structured explosion ejecta to be modelled using a static target in an implosion 
q | facility. 

1 Introduction 

There is much interest at present in the possible use of the new generation of high-power laser facilities 
(in particular the National Ignition Facility at Livermore and the Laser MegaJoule in Bordeaux) to 
simulate astrophysical phenomena such as supernovae. At first sight this programme appears to suffer 
from one obvious drawback. The phenomena one wishes to simulate generally involve explosions while 
the laser facilities are designed to produce implosions. Remarkably, as we will show, this is not a 
problem. Under certain, not too restrictive, conditions there exists an exact mathematical duality 
which allows one to transform an explosion problem to an implosion problem and vice versa. Thus it 
is possible, in a precise sense, to use implosion experiments to simulate exploding systems. 

2 The duality transformation 

The Euler equations of perfect gas dynamics can be conveniently written in the form 

-pV • U (1) 
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where D/Dt denotes the Lagrangian, material or convective derivative defined by 

and p is the mass density, £ the thermal energy density, p the pressure and U the velocity. In 
addition to these differential equations one algebraic relation is needed, an equation of state relating 
the pressure to the mass and energy densities the simplest being a polytropic equation of state, 

P= (1-1)8 (5) 

with 7 a constant. These equations are mathematically equivalent to the mass, momentum and energy 
conservation equations in smooth regions of the flow. Only at shocks is it necessary to revert to the 
fundamental conservation forms to recover the correct shock jump conditions. 

Now consider the following transformation of the dependent and independent variables, 

x* = a(t) _1 x, (6) 
t* = J a(t)~ 2 dt, (7) 

P* = a(tfp, (8) 

V = a(tfp, (9) 

IT = a(t)U-d(t)x, (10) 

£* = a{tf£ (11) 

where for the moment a(t) is an arbitrary function of time. Apart from the, at first sight rather 
strange, time-dependent scaling factors this is essentially a transformation to a coordinate system 
which is expanding or contracting with a scale factor a(t). If we define a* = a~ l and note that 

da* 2 d (1\ da 

~dt* ~ a ~dl (a) ~ ~~dt ^ ' 

it is easy to see that it is an involutionary transformation with inverse obtained by simply interchanging 
the starred and unstarred quantities 

x = a*(t*) _1 x*, (13) 
a*(t*)- 2 dt*, (14) 

P = a*(t*fp*, (15) 

p = a *(t*) 5 p*, (16) 

U = a*(t*)XJ* -d*(t*)x*, (17) 

£ = a*(t*f£*. (18) 

Let us now consider how the dynamical equations transform under this change of variables. It is 
easy to see that 

Dt Dt* v ; 

V -» a*V* (20) 
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and thus, after some elementary algebra, 



Dp* 
~Dt* 



-p*V*-U*, (21) 



mr = vy a* 

Dt* p* a* ' V ; 

(£* + p*)V* ■ U* + — (3p* -26*). (23) 



Remarkably, we see that if the scale factor a is such that 

... d 2 a* 



dt* 2 



a 2 d = (24) 



and the gas is a polytrope of exponent 5/3 with p = 28/3 then the Euler equations are invariant 
under this transformation. Note that, because the Euler equations in conservation form are alge- 
braically equivalent to the simplified forms, the conservation forms are also invariant and thus the 
whole structure of ideal gas dynamics, including the Rankine-Hugoniot shock relations, is preserved. 

The condition that the acceleration of the scale factor be zero, a = 0, requires that a(t) be a linear 
function of t and, without loss of generality, we can take a = t/to where to is a constant characteristic 
expansion time. The time transformation is then 



f dt o r dt tf. 



and it is convenient to set the constant to zero and choose 

(26) 



f 2 fl 
l . _ r 



The initial singularity of the expansion in physical space occurs at t = and is mapped to t* = — oo, 
the long term behaviour as t — > +oo is mapped to t* = 0. It is important to note that in the dual 
representation the time variable is bounded from above, t* < 0, whereas in physical space it is bounded 
from below, t > 0. 

The remarkable result is that for an ideal gas of point particles with no internal structure (which 
is what the 5/3 polytrope is) hydrodynamics in a uniformly expanding system is exactly equivalent to 
hydrodynamics in a static system. This result, or special forms of it, appears to have been discovered 
a number of times by Cosmologists (where the idea of factoring out the general expansion of the 
universe is very natural); a recent discussion is that of Martel and Shapiro (1998) where they propose 
the felicitous name of "supercomoving variables" to describe this transformation. What does not seem 
to have been generally noted is that this transformation can be used outside the cosmological context 
(however Poyet and Spiegel, 1979, did use a variant in an analysis of stellar pulsations). 



3 Interpretation 

The fact that the transformation is exact for the gas of ideal point particles strongly hints that it is 
derived from a similar result for the free particle motion. In fact there is such a duality, although it 
is almost trivial. The freely moving point particle moves along a straight line trajectory, 

x = x + v t, (27) 
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with starting point xq and velocity vq . If we write this as 



^ = Vo + Xo I (28) 

we see that there is a dual representation of the trajectory in which t is replaced by l/t, lengths are 
scaled by a factor proportional to time, initial points and final velocities are interchanged, but the 
trajectory remains a straight line. If collisions are instantaneous, localised and elastic they look the 
same in either system, and thus in both systems one can write down a Boltzmann equation and then 
derive the hydrodynamic equations as limits of moments of the Boltzmann equation. This approach 
also shows that higher order effects, such as viscosity and heat conduction, can formally be treated in 
the same way; however the resulting transformed transport coefficients will in general have unphysical 
time dependencies (for an application see Drury and Stewart, 1976). 

This analysis also shows that similar results will hold in different numbers of spatial dimensions, 
but the equation of state will have to correspond to the ideal gas in that number of dimensions. In d 
spatial dimensions it is easy to verify that the "super-comoving" transformation takes the form 

x* = a(t) -1 x, (29) 
t* = J a(t)' 2 dt, (30) 

P* = a(t) d p, (31) 

P* = a(t) d+2 p, (32) 

U* = a(t)U-d(t)x, (33) 

S* = a{t) d+2 8 (34) 

and that the Euler equations are invariant if the gas has a polytropic equation of state such that 

dp = 28 (35) 

corresponding to an adiabatic exponent 

7=1 + ^. (36) 

An interesting way of looking at this transformation (for which we are indebted to our colleague 
Etienne Parizot) is that it provides an analogue in spherical geometry to the freedom that Galilei 
transformations allow in planar geometry. If we are looking at a planar shock, it is often convenient to 
transform to a reference frame where the upstream medium, or the downstream medium, or the shock 
itself, appears stationary. In spherical systems one cannot apply Galilei boosts because the origin is 
fixed, however this transformation, by allowing one to take out an arbitrary uniform expansion, gives 
one much the same freedom. 



4 Application to a Supernova Remnant 

Computational studies of the evolution of a Supernova Remnant commonly start with initial conditions 
of dense pressure-free ejecta expanding ballistically away from the site of the explosion, which it is 
convenient to locate at the coordinate origin, and interacting with a stationary, or slowly moving, 
ambient medium of much lower density and negligible pressure. To illustrate the application of the 
duality transformation let us consider the simple, if somewhat artificial, case of uniform density ejecta 
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interacting with a uniform and stationary ambient medium in perfect spherical symmetry. Then the 
initial conditions correspond to 



p(r, t) 


= Po 


U(r,t) 


r 

V 


p(r, t) 


= 



(37) 

(38) 
(39) 

in the region r < Vot occupied by the ejecta (Vo is the maximum expansion speed of the ejecta) and 

p(r,t) = Po (40) 
U(r,t) = 0, (41) 
p(r,t) = (42) 

in the external (r 3> Vot) medium of constant density po- The sweep-up time tsw corresponds to the 
point where the ejecta, if expanding unimpeded, would have a density equal to the ambient medium. 

This defines the physical problem of expanding ejecta interacting with a stationary environment. 
Let us now consider the dual problem obtained by applying the transformation with scale factor 



a(t) = * (43) 
tsw 



Then the dependent variables transform as 



r* = tsw^, (44) 
t * = _§K (45) 

so that the explosion, which occurs at t = in physical problem, occurs at t* = — oo in the dual 
problem. Conversely the asymptotic evolution as t — > oo in the physical problem is mapped to the 
behaviour at t* = in the dual problem. 

The ejecta density in the dual problem is constant, 

p*(r*,t*) = a 3 p(r,t)= Po (46) 

and the velocity is zero, U* = 0, in r* < Votsw- However the ambient medium is now time-dependent 
with density, in the region r* S> Voisw 

t \ 3 f-t*^- 3 



P*(r*,t*) = (— p = — po (47) 
and velocity 

U*(r*,t*) = ^. (48) 

Thus in the dual problem we have stationary ejecta interacting with an imploding ambient medium 
whereas in the physical problem we have exploding ejecta interacting with a stationary ambient 
medium. Instead of the initial explosion at t = in the physical problem we have the final crunch at 
t* = in the dual problem. 

The evolution in physical space of the supernova remnant structure has been often discussed and is 
well-known (eg Truelove and McKee, 1999; Dwarkadas and Chevalier, 1998). At early times, t <C tsw> 
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the bulk of the ejecta expand ballistically except for a thin interaction region on the outside consisting 
of a forward shock running into the ambient medium, a zone of hot shocked ambient medium, a 
contact discontinuity, a zone of shocked ejecta and a reverse shock propagating slowly into the ejecta. 
At later times, when the mass of swept up ambient material becomes comparable to the ejecta mass, 
the reverse shock detaches itself from the contact discontinuity and implodes on the centre and the 
outer forward shock approximates the self-similar Sedov solution for a strong point explosion in a cold 
gas. 

In the dual system the interaction looks a little different, and in some ways is simpler. Initially 
we have the stationary sphere of high density material (which for convenience we continue to call 
the ejecta, although in the dual representation it has not been ejected but is simply sitting there) 
surrounded by a very low density converging flow. The inflowing gas has to decelerate at a shock 
which stands about 10% further out in radius than the edge of the ejecta. Writing for convenience 
r = —t*/tsw there is an exact similarity solution in which U* oc t _1 , p* oc r~ 3 and p* oc r -5 in the 
region external to the sphere of ejecta. This steeply rising pressure (oc r -5 ) drives the reverse shock 
into the ejecta and starts the implosion of the ejecta. 

At later times, as the ejecta collapse, the shock in the imploding ambient medium also moves 
inwards thereby reducing the rate of increase of the pressure. Transforming the Sedov solution to the 
dual system we see that the shock radius scales as 

r* oc r 3 / 5 (49) 

and the postshock pressure as 

P*oct~ 19 / 5 . (50) 

Figure 1 attempts to show schematically the relation between the two representations. We note in 
passing that the dual representation is also useful for analytic and numerical studies; this aspect will 
be explored in a companion paper (Dwarkadas and Drury, in preparation). 

5 Prospects for laboratory simulations 

The perfectly symmetric explosion is neither realistic nor especially interesting; it is the easiest case 
to analyse numerically and there is no reason to suppose that a laboratory simulation would yield 
any additional information. However reality is more complicated. It is clear that the ejecta emerging 
from real supernova explosions are highly nonuniform on a wide range of scales and that to calculate 
the resulting remnant evolution in three dimensions is likely to remain a computationally challenging 
problem for some considerable time (cf Arnett, 1999). 

The interesting implication of this work is that it should be possible with the new generation of 
implosion facilities to simulate precisely this problem, the interaction of highly structured ejecta with 
their surroundings including all the effects of spherical geometry. One can easily imagine constructing 
a solid target whose density distribution models the density distribution of the expanding ejecta. If 
this target is then used in an implosion experiment, and if the momentum loading on the surface is 
tailored to rise in the same manner as the pressure behind the forward shock in the dual system, a 
steep initial rise as r -5 decreasing to r~ 3 - 8 , the evolution of the internal structures including all the 
turbulent mixing, instabilities and shock formation, should be exactly replicated. 

We emphasise finally that the transformation discussed in this paper is additional to and com- 
plements the well-known linear scaling relations as excellently discussed by Ryutov et al (1999) in 
the astrophysical context, or Connor and Taylor (1977) in the plasma physics context. Dimensional 
similarity and scaling arguments are obviously central to any attempt at simulation on a laboratory 
scale of astrophysical systems, however precisely because they are very general and linear they cannot 
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t=0 

Figure 1: A schematic representation of the SNR structure and evolution as seen in physical space 
(left) and the dual space (right). The locations of the forward shock (dashed) , contact discontinuity 
(dotted) and reverse shock (dash-dotted) are indicated. 
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turn an explosion into an implosion. The remarkable nonlinear symmetry discussed in this paper is 
specific to the ideal gas equation of state, but subject to this constraint gives a powerful new degree of 
freedom in simulation studies by allowing an arbitrary uniform expansion or contraction to be factored 
out thereby transforming an explosion problem to an implosion one or vice versa. 
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